Local correlations and hole doping in NiO 
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Using a combination of ab initio bandstructure methods and dynamical mean-field theory we study 
the single-particle spectrum of the prototypical charge-transfer insulator NiO. Good agreement with 
photoemission and inverse-photoemission spectra is obtained for both stoichiometric and hole-doped 
systems. In spite of a large Ni-d spectral weight at the top of the valence band the doped holes 
are found to occupy mainly the ligand p orbitals. Moreover, high hole doping leads to a significant 
reconstruction of the single-particle spectrum accompanied by a filling of the correlation gap. 
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I. INTRODUCTION 

Already in 1937, at the outset of modern solid state 
physics, de Boer and Verwej^ drew attention to the sur- 
prising properties of materials with incompletely filled 
3d-bands, such as NiO. This observation prompted 
Peierls and Mott 2 - to discuss the interaction between the 
electrons. Ever since transition metal oxides (TMOs) 
were investigated intensively. This interest further in- 
creased when it was discovered that TMOs display an 
amazing multitude of ordering and electron correlation 
phenomena, including high temperature superconductiv- 
ity, colossal magnetoresistance and Mott metal-insulator 
transitions^. In the late 1950's MnO and NiO were taken 
as the textbook examples of antiferromagnets. However, 
when the importance of local Coulomb correlations in the 
transition metal d-shell was realized TMOs were consid- 
ered candidates for Mott insulators 2 -. In the mid 1980's 
Zaanen, Sawatzky and Allen (ZSA) introduced their clas- 
sification of TMOs and related compounds into Mott- 
Hubbard and charge-transfer (CT) systems 4 -. In the early 
TMOs the ligand p-band is located well below the transi- 
tion metal d-band and thus plays a minor role in the low 
energy dynamics. Such a case, called Mott-Hubbard sys- 
tem in the ZSA scheme, is well described by a multi-band 
Hubbard model. On the other hand, the late TMOs be- 
long to the CT type where the p-band is situated between 
the interaction split d-bands. A more general Hamilto- 
nian where the p-states are explicitly included is then 
needed, which can be viewed as a combination of multi- 
band Hubbard and Anderson lattice models. A major 
impulse for detailed investigations of CT systems, and 
especially of their hole doped regime, came with the dis- 
covery of high temperature superconductivity in cuprate 
perovskites. While the standard three-band Hamiltonian 
for cuprates 5 . contains only one d-orbital per lattice site, 
the description of cubic transition metal monoxides, the 
prominent member of which is NiO, requires the full set 
of d-orbitals. The latter are of interest not only for funda- 



mental research, but play an important role also in fields 
such as geophysics^. Furthermore, recent progress in high 
pressure experiments^ made the insulator-to-metal tran- 
sition in some TMOs accessible in the laboratory, provid- 
ing yet another stimulus for theoretical investigations. 

We report a computational study of NiO combining 
ab initio band structure calculations in the local den- 
sity approximation (LDA) with the dynamical mean-field 
theory (DMFT), an approach known as LDA+DMFT^. 
By treating the local correlations and the Ni 3d - O 2p 
hybridization on the same footing we provide a descrip- 
tion of the full valence and conduction band spectra of 
a CT system with strong hybridization. We will show 
that a good quantitative agreement with photoemission 
and inverse-photoemission data can be obtained thereby. 
This provides a solid foundation for the subsequent in- 
vestigation of hole doping, a question of broader interest 
mainly in the context of cuprates. It will be shown that 
the behavior of the doped holes clearly reveals the CT 
character of NiO. 

NiO is a type II anti-ferromagnet (TV = 523 K) with 
a magnetic moment of almost 2 /jlb and a large gap sur- 
viving well above TV- The standard LDA band theory 
predicts NiO to be a metal^, or an antiferromagnetic 
insulator— if spin polarization is allowed. A severe un- 
derestimation of the gap and the magnetic moment sug- 
gests, however, that the Slater antiferromagnetic state 
obtained within LDA does not describe the true nature 
of NiO. On the other hand exact diagonalization studies 
on small clusters were quite successful in describing the 
single- and two-particle spectra 1 ^, showing that the local 
Coulomb interactions are important. This made it clear 
that an explicit treatment of Coulomb interactions within 
the 3d shell is needed, and methods such as LDA+U 12 , 
self-interaction correction^, or GW 1 - were introduced. 
The static, orbitally dependent self-energy of LDA+U 
enforces a separation of the occupied and unoccupied ri- 
bands and thus opens a gap comparable to experiment. 
This in turn leads to a significant improvement of the 
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TABLE I: Orbital occupancies and the local moment on the 
Ni site for different hole dopings. 
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description of static properties such as the local moment 
or the lattice dynamics 15 . However, the LDA+U method 
is limited to an ordered state and does not yield the elec- 
tronic excitations and the effect of doping correctly. 

A systematic inclusion of dynamical correlations was 
made possible by dynamical mean-field theor y 16 i 17 ' 18 . 
Since its introduction DMFT proved to be a power- 
ful tool for the investigation of electronic systems with 
strong local correlations. In connection with band- 
structure methods, the LDA+DMFT scheme^ provides 
access to material specific single-particle spectra as well 
as more general correlation functions. Applications 
of LDA+DMFT so far were mostly limited to Mott- 
Hubbard systems, where the ligand states are integrated 
out before the correlation problem is solved. Recently, 
Ren et al^ applied this approach to NiO and were able 
to obtain a realistic gap and the near-gap spectra. How- 
ever, this approach takes into account only d-electrons, 
such that the orbital character of the valence and conduc- 
tion band are bound to be the same, the high frequency 
incoherent features in the valence band are missing, and 
the hole doping cannot be described realistically. In this 
work we go beyond such limitations by working with the 
ligand p states explicitly. 



II. COMPUTATIONAL DETAILS 

LDA+DMFT proceeds in two steps: (i) construction 
of the effective Hamiltonian from converged LDA cal- 
culation, and (ii) solution of the corresponding DMFT 
equations. Here we use the projection onto Wannier 
functions 2 ^ to obtain an eight-band p — d Hamiltonian 

h Z^ d i a aP^l" + K d -y a Pl-ya d ^a)+ (1) 

Here dkao- and pk 7tT are Fourier transforms of di aa and 
Pi-fa, which annihilate the d or p electron with orbital 
and spin indices ever or 7er in the ith unit cell, and nf aa is 
the corresponding occupation number operator. The ele- 
ments of matrix are parameterized by U and J. The 
constrained LDA calculation yields U=8 eV and J=l 
eV— . To account for the Coulomb interaction already 
present in LDA we renormalize the dd-diagonal elements 
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Experiment 




Energy (eV) 

FIG. 1: (color online) Theoretical Ni-d (solid line) and O-p 
(shaded) resolved spectral densities compared to photoemis- 
sion and inverse photoemission data obtained at 120 eV and 
66 eV photon frequencies after Ref. [21]. Gaussian broadening 
of 0.6 eV full width at half maximum corresponding to the 
experimental resolution was applied to the theoretical curves. 
The circle marks position of the d 10 L excitation. 

of the LDA Hamiltonian by the double counting correc- 
tion 

K d a p = ~K d a p{K) - (Novb - l)Un LDA 5 aP (2) 

where ti^da is the average LDA occupation per orbital 
and A r or b=10 is the total number of orbitals within the 
shell. 

Next we iteratively solve the DMFT equations on the 
Matsubara contour, a key part of which is the auxil- 
iary impurity problem solved by quantum Monte-Carlo 
(QMC) method 2 ^. The results reported here were ob- 
tained at T=1160 K. To obtain the single-particle spec- 
tral functions analytic continuation to real frequencies is 
performed using the maximum entropy method 23 . Ap- 
plying QMC to a gapped system requires careful assess- 
ment of ergodicity and autocorrelation issues 2 ^. 

III. RESULTS AND DISCUSSION 

A. Single particle spectra 

The orbital occupations, shown in Table [1] and the lo- 
cal moment of 1.85 fiB obtained from the paramagnetic 
DMFT solution correspond to a e? 8 groundstate of the Ni 
ion with two ferromagnetically coupled holes of e g sym- 
metry. In Fig. [T]the calculated spectral densities resolved 
into Ni 3d and O 2p contributions are compared to pho- 
toemission and inverse photoemission data3i. Using the 
full p — d Hamiltonian we are able to cover the entire 
valence and conduction bands spectra. Features corre- 
sponding to 4s and 4p bands at 10 eV and 13 eV, respec- 
tively, are not included in the theoretical spectrum. As 
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FIG. 2: (color online) Orbitally resolved O 2p and Ni 3d spec- 
tral densities (offset for better resolution). Substantial spec- 
tral weight transfer relative to the LDA results (see inset) is 
observed. 



FIG. 3: (color online) Ni-d and O-p resolved spectral densities 
for a hole concentration nj, = 0.6 (offset for better resolution). 
The inset shows a comparison of e g spectral densities for hole 
concentrations nn = 0.6 and 1.2. 



shown by Eastman and Freeouf 2 ^ the relative intensity 
of the 2p contribution increases with decreasing photon 
energy. Therefore the 120 eV spectrum is dominated by 
Ni 3d emission, while at 66 eV photon energy the O 2p 
contribution peaked around -4 eV is resolved (for a de- 
tailed orbital decomposition see Ref.—). The theoretical 
spectrum very well reproduces the experimental features, 
including the size of the gap, the d character of the con- 
duction band, the broad d peak at -9 eV, the position 
of the p-band, and the strong d contribution at the top 
of the valence band. While the gap and the Hubbard 
subbands can be described already with the static the- 
ory (LDA+U)^, a dynamical treatment is apparently 
needed to capture the substantial redistribution of spec- 
tral weight between the incoherent (-9 eV) and resonant 
(-2 eV) features in the d spectrum. For a more detailed 
analysis we show in Fig. [5] the spectral density resolved 
into e g and tig representations. In agreement with other 
studie d 11 ' 19 we find the conduction band to have a pure e g 
character. The dominant feature of the valence spectrum 
is a distribution of spectral weight between the broad 
peak at high frequency and sharp peak(s) at the gap 
edge. The origin of these structures was discussed by Fu- 
jimori et aZ»ii in terms of eigenstates of a NiC>6 cluster. 
Emitting a d electron from the d s state the system can 
end up either directly in a d 7 final state (broad peak) 
or in a d a L final state (sharp peak), with a ligand hole, 
due to p — d electron transfer. We use this picture to 
discuss our DMFT results. A sizable d 8 L_ resonant peak 
appears in the filled tig band, in spite of only a weak 
V ~ ^2g hybridization. This is interpreted as being due 
to the emission of a tig electron followed by an electron 
transfer from oxygen to the partially filled e g shell. The 
e g spectrum exhibits even richer structure. In the con- 
duction sector we observe, besides the dominant d 9 peak, 
also a tiny high frequency feature due to the d 10 L final 



state, an accurate resolution of which is difficult using 
the maximum entropy method. The the resonant d s L 
peak exhibits a pronounced splitting in the e g channel. 
We ascribe this to formation of the Zhang-Rice bound 
stated at the gap edge, which was studied in detail in 
Ref. Similar effect was also observed on a simpler 
two band p-d mode l 28 ' 29 . In the present work an accu- 
rate resolution of the fine structure of the resonant peak 
is not possible since the spin-flip terms in the interaction 
are not included. 



B. Hole doping 

Next we discuss hole doping of NiO. An experimen- 
tal realization can be found in Li^Nii-zO studied in the 
doping range a;=0.02-0.4 30 . Using the Hamiltonian of 
the undoped system the replacement of x Ni 2+ ions by 
Li 1+ ions introduces on average rih = x/(l — x) hole per 
Ni site. As crude as this approximation may be we be- 
lieve that the essential physics of p-d weight transfer is 
captured correctly. In Fig. [3] we show the single-particle 
spectral densities for = 0.6 corresponding to x = 0.38. 
There is no significant difference between the ti g spectra 
in the doped and undoped cases, but the e g spectral func- 
tion changes significantly. Most notably the Mott gap is 
filled, while the Hubbard subbands are preserved as dis- 
tinct features. This is also observed in experiment, as 
shown in Fig. [4] A quite different behavior was reported 
in the single-band Mott insulator—, where the gap sur- 
vives doping while the coherent peak merges with one 
of the Hubbard subbands. Further hole doping of NiO 
leads to a spectral weight transfer from both upper- and 
lower-Hubbard-subbands to the quasiparticle part of the 
e g spectrum (see inset of Fig. [3]), which can be viewed 
as an enhancement the itinerant character of the system 
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FIG. 4: (color online) Theoretical Ni-d spectral densities for 
electron addition and electron removal obtained for nn = 
0.6 hole-doped NiO compared to the photoemission and in- 
verse photoemission spectral of Lio.4Nio.6O (the experimen- 
tal baseline is offset for better readability). 

reflected also in the decrease of the local moment (see Ta- 
ble UJ). The p spectral density suggests that much of the 
doped hole is accommodated on the oxygen sites. That 
this is indeed the case can be seen from the reduction 
of the p orbital occupancy shown in Table |TJ In an un- 
corrected system doping results in a mere shift of the 
chemical potential and the variation of the orbital occu- 
pancy is given by the spectral density at the chemical 
potential. This reasoning clearly fails in a multi-orbital 
correlated system such as NiO. While the Ni-cZ spectral 
density at the top of the valence band in the stoichio- 
metric system (or at the chemical potential in the doped 
system) is comparable to the O-p spectral density, doped 



holes almost exclusively reside at the oxygen sites. 



IV. CONCLUSIONS 

We used the LDA+DMFT(QMC) approach to com- 
pute single-particle spectra of the prototypical CT in- 
sulator NiO. By including the ligand p states and the 
on-site Coulomb interaction within the same framework 
we were able provide a full description of the valence 
band spectrum and, in particular, of the distribution 
of spectral weight between the lower Hubbard band 
and the resonant peak at the top of the valence band. 
Good agreement with the available photoemission and 
inverse-photoemission data was found without need for 
adjustable parameters. Importantly, the present method 
allows us to study the hole doped regime where we find 
that the doped holes are distributed mainly among the 
ligand sites. A high hole doping leads to the filling of the 
correlation gap and a significant transfer of the <i-spectral 
weight form the incoherent part of the spectrum. 
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